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Abstract. We present the first evolutionary models of 
intermediate mass stars up to their thermal pulses which 
include effects of rotation on the stellar structure as well as 
rotationally induced mixing of chemical species and angu- 
lar momentum. We find a significant angular momentum 
transport from the core to the hydrogen-rich envelope and 
obtain a white dwarf rotation rate comparable to current 
observational upper limits of < 50 kms^^. 

Large angular momentum gradients at the bottom of 
the convective envelope and the tip of the pulse driven 
convective shell are shown to produce marked chemical 
mixing between the proton- rich and the ^^C-rich layers 
during the so called third dredge-up. This leads to a sub- 
sequent production of ^^C which is followed by neutron 
production through ^^C(a, n) in radiative layers in be- 
tween thermal pulses. Although uncertainties in the effi- 
ciency of rotational mixing processes persist, we conclude 
that rotation is capable of producing a ^^C-rich layer as 
required for the occurrence of the s-process in TP-AGB 
stars. 

Key words: Stars: evolution - Stars: abundances - Stars: 
AGB, post-AGB - Stars: rotation 



1. Introduction 

It is known since long that thermally pulsing Asymp- 
totic Giant Branch (TP-AGB) stars provide a site for 
the so called s-process, i.e., the slow neutron capture 
process which forms neutron-rich isotopes heavier than 
iron (Clayton 1968). Heavy elements primarily produced 
by the s-process are overabundant at the surface of 
AGB stars (Smith & Lambert 1990), including tech- 
netium (Little et al. 1987) which has no stable isotope 
and which is produced as ^^Tc (ti/2 = 2.1 x 10^ yr) 
in the s-process. In particular the roughly solar magne- 
sium isotopic pattern found in s-process enriched AGB 
stars has demonstrated that the "'^'^C(a,n) rather than the 
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^^Nc(a,n) neutron source is likely to operate the s-process 
in AGB stars (Guehn et al. 1995, Lambert et al. 1995). 

Evidence for in situ s-processing is found exclusively 
in carbon stars (Smith & Lambert 1990), which corre- 
spond to a late evolutionary stage on the TP-AGB where 
the stars have large ^^C enrichments in their envelopes 
(Iben & Renzini 1983, Wallerstein & Knapp 1998). The 
^^C enrichment implies that these stars contain, at cer- 
tain times, a region at the bottom of their hydrogen- 
rich envelope where ^^C is abundant. This region where 
protons and ^^C coexist may then perhaps form ^■^C 
through i2C(p,7)^^N(/3+i/)i3C. Although this scenario is 
unrivaled, the formation of a layer which is rich in protons 
and ^^C in TP-AGB models has proven to be difficult, and 
its existence had hitherto to be assumed ad hoc in all s- 
process calculations (Gallino et al. 1998). Iben & Renzini 
(1982) found a -'^'^C layer in low metalicity ABG models. 
Recently, Herwig et al. (1997) have obtained a ^'^C-rich 
layer in TP-AGB models of solar metallicity, by invoking 
a diffusive overshoot layer at convective boundaries. Here, 
we investigate for the first time effects of rotationally in- 
duced mixing processes on the TP-AGB. 

2. Numerical method and physical assumptions 

Our calculations have been performed with a hydrody- 
namic stellar evolution code (cf., Langer 1998, and refer- 
ences therein), which has been upgraded to include angu- 
lar momentum, the effect of the centrifugal force on the 
stellar structure, and rotationally induced transport of an- 
gular momentum and chemical species due to Eddington- 
Sweet circulations, the Solberg-H0iland and Goldreich- 
Schubert-Fricke instability, and the dynamical and secular 
shear instability. We apply the rotational physics exactly 
as in Heger et al. (1999). In particular, the effects of gradi- 
ents of the mean molecular weight /x, which pose barriers 
to any mixing process, have been included as in Heger 
et al. (i.e., = 0.05). As in Heger et al., we have also 
included the effects of /^-barriers on convection by using 
the Ledoux-criterion for convection and semiconvection 
according to Langer et al. (1983), which is consistent with 
our treatment of the rotational mixing (Maeder 1997). 
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Fig. 1. Evolutionary track of our rotating 3Mq model 
(solid line) and of a non-rotating reference model (dotted 
line) in the HR diagram. The tracks start at the zero age 
main sequence and end on the TP-AGB. The beginning 
of the thermal pulses is marked. 

Changes of the chemical composition and the nuclear 
energy generation rate are computed using nuclear net- 
works for the three pp-chains, the four CNO-cycles, and 
the NcNa- and the MgAl-cyclc. Further, the 3a-reaction is 
included, and (Q;,7)-reactions on ^^C, ".isn^ le.iSQ, Wp^ 
20,21,22^6, 24,25,26Mg, and (a,n)-reactions on ^^C, ^^O, 
2i,22^g^ 25,26]y[g^ rpj^^ inclusion of (n,7)-rcactions on -^^C, 
^°'^^Ne, ^"^'^^Mg, 28,29g- allows an order of magnitude es- 
timate of the neutron concentration. For more details see 
Heger et al. (1999) and Heger (1998). 

3. Results 

3.1. Evolution towards the TP-AGB 

For our pilot study, we chose to compute the evolution of 
a 3Mq star of roughly solar composition with an initial 
equatorial rotation velocity of 250kms~^, which is typical 
for late B main sequence stars (Fukuda 1982). This choice 
renders effects of magnetic braking and the core helium 
flash unimportant. The evolution of our model in the HR 
diagram, together with that of a non-rotating reference 
model, is shown in Fig. 1. 

As in massive stars (Heger ct al. 1999), the dominant 
rotational mixing process on the main sequence is the 
Eddington-Sweet circulation. It leads to a ^^C/^^C-ratio 
after the first dredge-up of 9.4, compared to a value of 
19.4 in our non-rotating model (cf., Boothroyd & Sack- 
mann 1999). 

Figure 2 sketches the evolution of the angular momen- 
tum distribution in the innermost 0.8 Mq of our rotating 




I 1 1 1 I I I i I , i_L_i I L. I I I 

0.2 0.4 0.6 0.8 

K / Mo 

Fig. 2. Upper panel: logarithm of the local specific an- 
gular momentum (in cm^ s~^) as function of the mass co- 
ordinate for three models of our rotating 3 Mq sequence; 
at core hydrogen exhaustion {t ~ 3.3 x 10^ yr, solid line), 
during core helium burning {t ~ 3.8 x 10^ yr, dotted 
line), and between the 14th and the 15th thermal pulse 
{t ~ 4.4 X 10^ yr, dashed line). Lower pEinel: logarithm 
of the local angular velocity (in rad/s) as function of the 
mass coordinate for the same models which are displayed 
in the upper panel. 



3M0 model. It shows that the core specific angular mo- 
mentum decreases continuously during the evolution. We 
can give a first quantitative prediction of a white dwarf 
rotation rate: log(77cm^s^^) = 15.3 and i?,wD = 0.01 R© 
yields fj-ot = 28kms^^. This is of the same order as cur- 
rent observational upper limits (Heber et al. 1997, Koester 
et al. 1998). 

3.2. Mixing and nucleosynthesis on the TP-AGB 

Figure 3 shows the evolution of the internal structure dur- 
ing and after the 25th thermal pulse of our rotating model. 

It shows that the tip of the pulse-driven convection zone 
leaves after its decay a region of strong rotational mix- 
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Fig. 3. Section of tlie internal structure during and af- 
ter the 25th thermal pulse of our rotating 3 Mq sequence. 
Diagonal hatching denotes convection. The convective en- 
velope extends down to Mr ~ 0.746 Mq. The pulse driven 
convective shell is located at 0.737 M© < < 0.746 M© 
and 30yr < t < 120 yr. Vertical hatching denotes regions 
of significant nuclear energy generation, i.e., the hydrogen 
burning shell (at Mr ~ 0.746 and t < 100 yr) and the 
helium burning sheU (0.734 Mq < Mr < 0.739 M© and 
t > 40 yr). Gray shading marks regions of significant rota- 
tionally induced mixing (see scale on the right side of the 
figure). Vertical marks at the bottom of the figure denote 
the time resolution of the calculation, where every fifth 
time step is indicated. Cf. also Fig. 4. During this thermal 
pulse, the maximum energy generation rate of the helium 
burning shell was 4 x 1O''L0. 

ing. This mixing becomes even stronger when the convec- 
tive envelope extends downward during the third dredge- 
up event (cf. also Figure 4). The reason is that convec- 
tion enforces close-to-rigid rotation (cf. Heger et al. 1999), 
with an envelope rotation rate which is many orders of 
magnitude smaller than that of the core. The resulting 
strongly differential rotation (cf. also Figure 2) allows the 
Goldreich-Schubert-Fricke instability, and to a lesser ex- 
tent the shear instability and Eddington-Sweet circula- 
tions, to produce a considerable amount of mixing between 
the carbon-rich layer and the hydrogen envelope. 

Figure 5 depicts the resulting hydrogen and ^^C abun- 
dance profiles after the convective envelope has receded. 
It shows a layer of several 10~^ M© containing a large 
mass fraction of protons and ^^C at the same time. Sev- 
eral 1000 yr after the pulse, this layer heats up and ^'^C is 
formed through proton capture on ^^C. Figure 5 shows the 
resulting ^"^C profiles for four different times. A maximum 
^^C mass fraction of almost 4% is achieved. 

Starting some 10^ yr after the pulse, the ^"^C-rich 
layer becomes hot enough for a-captures on ^^C to oc- 
cur (Straniero et al. 1995). Figure 6 shows the resulting 
neutron densities for three different times. Note that 
we did not include the reaction ^''N(n,p)^''C in our net- 
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Fig. 4. Same as Figure 3, for the same time interval, but 

magnifying the dredge-up of the convective envelope. Note 
the hydrogen burning shell source at Mr ~ 0.7461 and t < 
100 yr, and the extension of the pulse driven convection 
zone up to Mr ~ 0.74602 at t ~ 100 yr. 

work. Although most of the resulting protons may form 
new ^^C, it may be an effective neutron sink (Jorissen & 
Arnould 1989), in particular as also a large abundance 
of (primary) ^"'N is produced in the ^^C-rich layer. Thus, 
our neutron densities can only be considered as an order 
of magnitude estimate. With rin — 10^ cm^'"* for ~ 10"* yr, 
we obtain a neutron irradiation of r ~ 10^^ neutrons/cm^ 
which results roughly in a number of neutron captures per 
iron seed of ric ~ 75, i.e. a main component s-process (cf. 
Figures 7.22 and 7.23 of Clayton, 1968). 

4. Discussion 

By applying the concept of rotationally induced mixing as 
it has been developed for massive stars in our group dur- 
ing the last years without alteration to a 3 Mq TP-AGB 
model sequence, we obtain conditions which appear favor- 
able for the development of the s-process, i.e. a ^^C-rich 
layer which produces a considerable neutron flux later-on. 
Although our model develops only a very late and weak 
third dredge-up we believe that the mechanism which dif- 
fuses the protons into the carbon layer and -"^^C into the 
envelope must occur with a similar magnitude in all TP- 
AGB stars which develop a third dredge-up. The reason is 
that the huge specific angular momentum jump at the hy- 
drogen/carbon interface — five orders of magnitude in our 
case — is independent of the depth of the third dredge-up. 

The maximum "'^^C abundance and its distribution in 
our model is, at first, similar to that found due to dif- 
fusive convective overshooting by Herwig et al. (1997). 
However, in our case the rotational mixing spreads the 
^^C peak out before the neutrons are produced (cf. Fig- 
ures 5 and 6), which is not the case in the models of Herwig 
et al. (1997). At the present time we can not discriminate 
which of these scenarios would agree better with empirical 
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Fig. 5. Chemical profiles at the location of the maximum 
depth of the convective envelope during the 25th thermal 
pulse (cf. Figure 4) of our rotating SM© sequence. The 
dotted and dashed lines mark the hydrogen and the ^^C 
mass fractions at t = 1704 yr, with t — ^ defined as in 
Figures 3 and 4. The fat solid line denotes the ^^C mass 
fraction at the same time. The three thin solid lines repre- 
sent the ^^C mass fractions at t = 2016 yr, t = 4155 yr, and 
t = 5139 yr, with a later time corresponding to a smaller 
peak abundance. The maximum ^^C mass fractions of 
3.6% occurs at t = 2016 yr. The ^^C peak moved inwards 
in the time interval from t = 1704 yr to t — 2016 yr due to 
continued proton captures on both, ^^C and ^■^C. 



constraints. However, we want to stress that both mech- 
anisms of ^^C production, rotation and overshooting, do 
not exclude each other, and that it is possible that they 
act simultaneously in AGB stars. 

Finally, we want to emphasize that, although stars of 
less than ~ 1.3 Mq lose 99% of their angular momentum 
due to a magnetic wind during their main sequence evolu- 
tion, it can not be excluded that the proposed mechanism 
of ^'^C-production due to differential rotation also works 
for them. Certainly, the sun's core will spin-up and the 
envelope will further spin down during its post-main se- 
quence evolution, which may result in a specific angular 
momentum jump of similar magnitude. The investigation 
of the mass and metallicity dependence of the production 
of ^^C due to rotation is an exciting task for the near 
future. 
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Fig. 6. Neutron density as function of the mass coordinate 
for three models of our rotating 3 Mq sequence after the 
25th pulse. The fat solid line corresponds to t = 14 150 yr, 
the other two (thin solid lines) to t = 11457yr (peak at 
largest mass coordinate) and t — 16 727yr, with t — 
defined as in Figures 3, 4 and 5. The time span between 
the 25th and the 26th thermal pulse is ^ 30 000yr. 
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